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Abstract 

This paper concerns quantitative analysis of errors generated by in¬ 
completely known data in convex minimization problem^J. The problems 
are discussed in the mixed setting and the duality gap is used as the 
fundamental error measure. The influence of the indeterminate data is 
measured using the worst case scenario approach. The worst case error is 
decomposed into two computable quantities, which allows the quantitative 
comparison between errors resulting from the inaccuracy of the approxi¬ 
mation and the data uncertainty. The proposed approach is demonstrated 
on a paradigm of a nonlinear reaction-diffusion problem together with nu¬ 
merical examples. 


1 Introduction 

In mathematical models, the parameters (data) are often not completely known. 
Instead, their values are based on measurements or assumptions. For analysts, 
it is important to quantify the effect that the data uncertainty has on results 
generated by the mathematical model. In particular, the data uncertainty gen¬ 
erates an accuracy limit, beyond which information cannot be extracted and all 
(numerical or analytical) efforts become meaningless or even misleading. 

The analysis of uncertain data in PDFs (in the context discussed in this pa¬ 
per) dates back to [2] (see also and references therein). One of the main 

lines of the uncertainty analysis is the probabilistic approach, where the uncer¬ 
tain data is treated as a random variable. This leads to stochastic PDFs (see, 
eg., lanni and references therein). A recent overview of the related numerical 
methods can be found from [TH]. Analysis of stochastic PDFs with additional 
noise can be found from HMD- 

The analysis in this paper is based on the duality theory (convex analysis) 
developed by W. Fenchel, J.-J. Moreau, and R. T. Rockafellar (see, e.g., [inii^ 
and references therein). S. Repin applied the duality theory to derive com¬ 
putable a posteriori error estimates for convex problems in ^511251125] . In this 
paper, a class of convex variational problems (cf. (12.11) 1 are discussed and the 
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quantity of interest is the maximal duality gap related to a family of problems 
generated by the admissible (incompletely known) data. The approach is mo¬ 
tivated by the worst case scenario method (see, e.g., M and [21 Ch. 6.4]). In 
contrast to the probabilistic methods, the worst case scenario method does not 
assume which data are more or less likely than other, but is solely interested 
to find the most drastic (measured in a problem dependent manner) event pos¬ 
sible. The benefit of the worst case scenario approach is that there is no need 
to make any assumptions concerning the probability distribution of unknown 
parameters, but the results may be overpessimistic. 

In this paper, computable quantities (or two-sided estimates) related to the 
approximation error and the error generated by the incompletely known data 
are presented. This is of practical importance. Firstly, it is meaningless to use 
computational resources to improve approximations which are “too close” to (or 
even inside) the solution set. Secondly, it it crucial to evaluate whether the 
indeterminacy of the data is small enough to justify any implications based on 
the results (regardless of any numerical aspects). 

For linear problems, the main ideas for the uncertainty analysis (related to 
the one presented in this paper) were first discussed in [231 Ch. 6.8] and [271 Ch. 
9.6]. These ideas rely on the use of functional a posteriori error estimates, which 
are guaranteed and depend explicitly on the problem data. Additional research 
is exposed in [Mia. It is worth mentioning that the results partially coincide 
with the results of [3] obtained in a probabilistic framework (see [301 Rem. 3.1]). 

The paper is composed as follows. In Sect. [H standard definitions and 
results from the duality theory are promptly exposed. The relevant part for 
the later analysis are, conjugate functionals, the definitions of the primal and 
dual problems, the duality gap, and the necessary conditions (subdifferential 
inclusions) in terms of the compound functionals (see, e.g., [l^). The use of the 
duality gap as the fundamental error quantity is motivated in Sect. 12.21 where 
it is shown that in linear examples (generated by quadratic energy functionals) 
certain important a posteriori error estimation results, e.g., the Prager-Synge 
estimate m, are representations of the duality gap. 

Analysis of the respective mathematical models, where the data are not com¬ 
pletely known is discussed in Sect. [3| The general ideas based on the duality 
theory are explained in Sect. 13.11 In Sect. 13.21 the new error quantification 
method is applied to a nonlinear reaction diffusion problem. Two-sided bounds 
for the error generated by the incompletely known data are presented in Theo- 
rem l3.II The practicality of these bounds is demonstrated in numerical examples 
in Sect. 13.31 


2 Convex variational problems 

2.1 Definitions 

Let V and V be a Banach spaces with V* and Y* as the corresponding dual 
spaces. The duality pairing is denoted by brackets and subindices, e.g., {v*,v)v, 
where v* G V* and v G V. The set of functions V K, which are pointwise 
supremum over a set of affine functions is denoted by r(y). Fo(R) is the 
subset of functions other than constants -l-oo and —oo. Henceforth, assume 
that F G ro(R) and G G Fq)!"). 
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The problems discussed in this paper are of the following type. The primal 
problem is 

J(m) = inf J{v) := F{v) + G{Av), (2.1) 

vSiV 

and the respective dual problem is 

I*{p*) = sup r{y*) := -F*{A*y*) - G*(-y*), (2.2) 

where F* and G* are the respective Fenchel conjugates 

F*{v*) := sup{{v*,v)v - F{v)} and G*(y*) := sup{{y*,y)Y - G{y)}, 

v&V vGY 

and A € C{V, Y) is a bounded linear operator with an adjoint A* G C{Y*, V*), 
i.e., 

{y*,Av)Y = {A*y*,v)v, Vu G F, y* G FA (2.3) 

The solutions of the primal and dual problem are characterized by the following 
Theorem (see, e.g., pm Rem 4.2] or [5^1 Thr 2.2]). 

Theorem 2.1. Let J be coercive on V. If there exists uq G V such that 
F{uo) < + 00 , G(Auo) < + 00 , and G is continuous at Auq, then there exists 
u GV and p* G Y* , which are solutions of problems S2.1\) and S2.S\) . respec¬ 
tively. Moreover, the duality gap (denoted by E(u,p*)) vanishes, i.e., 


E{u,p*):=J{u)-F{p*)=Q, (2.4) 

and 

A*p* G dF{u), (2.5) 

-p* G dG{Au). (2.6) 

In the problems (ED and (1^ . the most fundamental quantities are the 

primal and dual energies (J and I*). If an approximation {v,y*) gV xY* 
is inaccurate in the sense of these energies, then it is clear that drawing any 
conclusions based on it is highly dubious. The approximation error is naturally 
measured by the value of the respective duality gap E{v,y*). 

The duality gap can be studied in greater detail based on how well the 
approximation satisfies the subgradient conditions ()2.5|l and (12.611 . The following 
Proposition (see, e.g., m Prop. 5.1]) characterizes the set of subgradients. 

Proposition 2.1. Under the earlier assumptions, v* G dE{v) if and only if 

Df{v,v*) := F{v)FF*{v*) - {v*,v)v = 0. 

The functional Dp : V x V* ^ M. is called the compound functional. The 
Fenchel-Young inequality 

Df{v,v*)>0, 'IvGV,v*gV* 

and ProDOsition l2.1l show that Df{v, v*) is a reasonable measure of the violation 
of the relation v* G dF{v) ([211 Thr. 2.2] and [221 (7.1.3)]). 
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Remark 2.1. The duality gap Ii2.4\l can be decomposed as follows: 


Eiv, y*) = F{y) + G{Kv) + F*{A*y*) + G*{-y*) 

= Fiv) + G{Av) + F*{A*y*) + G*{-y*) - {-y\Av)Y - {A*y\v)v 

= DF{v,A*y*) + DG{Av,-y*), (2.7) 

where the two compound functionals measure violations of relations (dS) and 
dS), respectively. For analyst, the form in d3( is more practical than 
since the values of the primal and dual energies attained at the exact solution 
are unknown, but the compound functional values are always positive and vanish 
only at the exact solution. 

2.2 Examples 

2.2.1 Linear elliptic problem 

Let V and Y be Hilbert spaces endowed with the inner products (■:-W and 
and norms || • ||v and || • ||j/, respectively. Operators A : Y ^ Y and 
B : V ^ V are linear, bounded, symmetric, and strictly positive definite. Let 
£ G V*. This example concerns the linear problem 

{A*AA + B)u = £. (2.8) 

The primal energy functional associated with (12.8|) is 

Jiv) = ^\\Avr^ + ^\\v\\l-{£,v)v, 

where 

II 2 /IU := Vi-^y^y)Y and ||u||£5 := \/ {Bv,v)v 

are equivalent norms in Y and V, respectively. To present J(v) in the form 
(12.11) . we set 

F{v) := ^\\v\\l-{£,v)v and G{y) := ^\\y\\% 

The corresponding conjugate functionals are 

F*{v*) = h\\v*+^\\l^^ and G* (y*) = WWa-y 

where 

II 2 /IU -1 := \/{-^~^y,y)Y and ||'y||r5-i := \/{B-^v,v)v. 

Remark 2.2. Formally, one should write 

F*{v*) = k\\^y\v*+£)\\l^, and G^iy*) = 

where Ty and Xy denote canonical isomorphisms (mapping of the Riesz repre¬ 
sentation theorem for Flilbert spaces) V ^ V* and Y -A Y*. Hereafter, we omit 
this detail to lighten the notation. 

The compound functionals associated with relations (1^ and (1^ are 

DF{v,A*y*) = ^\\v\\l - {£,v)y + i||AV - {v,^*y*)v 

= l\\Bv-£-A*y*\\l^, (2.9) 
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and 


Dg{Av, -y*) = ^WKvWl, + -y*)y 

= \\\AKv + y*\\l,-r, ( 2 . 10 ) 

respectively. The fact that the solution {u,p*) satishes Df{u, A*p*) = 0 and 
Dg{Au, —p*) = 0 leads to the decomposed form of (|2.8II 

-Ap* + Bu = e, (2.11) 

p = -AAu. (2.12) 

Remark 2.3. The compound functionals are squared (weighted) L^-norms of 
the residuals related to the mixed form of the problem. Minimizing the duality 
gap w.r.t. primal and dual variables over some finite dimensional suhspaces 
leads to, e.g., least square finite element methods (see, e.g., and references 
therein). 

Remark 2.4. Substituting conditions \2.11\) and {2.1 ‘A) to Df(v, A*y*) and 
Dg{v, A*y*), respectively, yields 

DF{v,A*y*) = i||S(u - u) + A*{p* - y*)|||_, 

and 

DG{Av,-y*) = \\\AA{v-u) + y* -p*\W-i. 

Expanding the terms and summing provides an alternative representation for 
the duality gap (cross-terms cancel each other due to {2.3\) and {2.1 A) ) 


£;(u,y*) = i||u-u||2+i||A*(p*-j/*)||2„, + i||A(n-u)||^ + i||2/*-p1lVi. 

By {2. Til . {2.91) . {2.101) . and {2.12\) one obtains the so called error equality 


||u - uWl + ||A*(p* - y*)\\l-^ + ||A(u - u)\\\ + \\y*-p*\\\-. 

= ||S^_£_AVI||-i + MAz; + y*||Vi, 

which has importance in a posteriori error control and adaptive finite element 

methods (see mm)- 

U B = 0, then (12.8|) reduces to 

A*AAu = £. (2.13) 


The primal energy functional is 

J{v) = i||Au||3t - 


where 

F{v) :=-(e,v)v and G{y) := ^\\y\\l^. 
The corresponding conjugate functionals are 


F*{v*) 


0 , 

+ 00 , 


ifv*+i = 0, 

else, 


and G*(y*):=i||y*||Vi- 


5 












The respective compound functionals are 



0, iiA*y*+i = 0, 

+ 00 , else 


and 



which clearly correspond to the equilibrium condition A*p* + £ = 0 and the 
duality condition p* + AAu = 0, respectively. 

Remark 2.5. Assume that y* is equilibrated, i.e., A*y* + £ = 0. Then, the 


duality gap is simply Dq{Av, —y*). On the other hand, the duality condition 
for (u,p*) provides (again, by i2.3\) ) 

||.4Aw - y*\\\-i = ||-4Ad - AAu + p* - y*\\'\-i 


= ||A(u - uWa + \\P* - 2/11^1 + 2(A(z; - u),p* - y*)y, 


where the last term vanishes, since y* is equilibrated. This leads to the well 
known Prager-Synge estimate |24| (for equilibrated fluxes) 



2.2.2 The power growth functional with lower order terms 

Let V = V* = Y = LP(fl,R^), V* = Lflfl,R^), 

A = V, and A* = —div. Let p' > 1 and p > 1, q' and q are the conjugate 
numbers defined by ^ ^ = 1 and ^ ^ = !> respectively. The domain 

n C R'^ is open, bounded and has a piecewise Lipschitz continuous boundary. 
Moreover, let p < p', then (fl,R) C LP{n,R) (see, e.g., |31l Th. 2]) and 
Wq '^ (n, K) C LP{n, R). This example concerns a p-Poisson problem with lower 
order terms, i.e.. 



M = 0, on did. 


These Euler’s equations are related to the primal energy functional 



(2.14) 


n 


where / S Lfln,R) and 7,7' e L^{n,R+). Here, 



and 



n 


The corresponding conjugates are 



and 



n 


n 
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and the compound functionals are (assume henceforth additional regularity 
divy* e 

Df{v, -divy*) = J ^ ( 7 ) I “ divy* + /|« - (-divy* + f)v^ dx 

a 

and 

DG{Vv-y*)= j (^^y|Vur' + i(i)'”'|y*r'-Vr;-(-y*)) dx. 

Q 

Remark 2.6. Note that in the case p = p' = 2 (then also q = q' = 2) 

Df{v, -divy*) = \ J ^ ( 71 ; + divy* - ff dx 
a 

and 

DoiVv, -y*) = lj± ( 7 'Vu + y*)' dx. 
n 

Assuming enough regularity, conditions Df{u, —divp*) = 0 and Dg{Nu, —p*) = 0 
lead to 


divp* + yu = /, in O, 

p* = — 7 'Vm, in Vt, 
u = 0, on d^l. 

3 Incompletely known data 

3.1 Worst case scenario approach 

In many mathematical models functionals F and/or G depend on some incom¬ 
pletely known data or parameters. Henceforth the set of these parameters is 
denoted by T. Motivated by the engineering practise, we assume that the set 
of admissible data is of the form 

T := {7 = 7 o -h(5, ||(5 ||t < e}, 

where 70 are the known “mean data” and 5 are the respective (bounded by some 
suitable norm || • ||x) variations resulting, e.g., from the measurement inaccuracy, 
incomplete known material parameters, etc. It is assumed that problems (EH) 
and (EH) have solutions for all admissible data. This relation is denoted by 

<S:T^<S(T) cHxF*. 

Note that every 7 G T generates not only a problem and the associated solu¬ 
tion 5 ( 7 ), but also the proper error measure (i.e., the duality gap) E-y (subscript 
denotes the dependence of functionals on the parameter 7 ) associated with that 
problem. The fact that the error measures are not uniquely defined stimu¬ 
lates the following definitions of extremal errors associated with an arbitrary 
{v,y*)&VxY\ 

E^°^^\v,y*)-.= sxxpE^[v,y*) 

7GT 
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and 


E^-^\v,y*) := ini E^{v,y*). 

7GT 

Remark 3.1. Recall that by i2. 71) . the duality gap is a sum of compound func¬ 
tionals 

E-f{v,y*) = DF^{v,A*y*) + Dg^{Ev,- y*). 

Since compound functionals vanish at the exact solution (and are positive oth¬ 
erwise), they could have been used independently to define similar worst (and 
best) case errors. 

Due to the computational cost, typically simulations are performed using 
only the mean data values. Thus the objective of computations is to approxi¬ 
mate the mean solution {uo,p*) '■= <S( 7 o). In the context of the mixed setting 
discussed in this paper, the approximation is “good” if Ej^ {v, y*) is small. (Here¬ 
after, we denote E^^ by Eo to avoid double subscripts.) However, if the data 
is not completely known, then there is an accuracy limit beyond which further 
computations (e.g., adaptive mesh refinements) make no sense. The following 
analysis provides tools to compare between the error resulting from the inaccu¬ 
racy of the approximation and the error generated by the data uncertainty. 

Assume that it is possible to explicitly linearize the duality gap functional 
with respect to the data perturbation 


Ej{v,y*) = Eo{v,y*) -f R{v,y*]-fo,S), 


(3.1) 


where R is the remainder term related to the particular mean data 70 and 
perturbation 6. Then can be decomposed by (EU into two parts 

= E,iv,y*)+E^’^^iv,y*), (3.2) 

where 

£''"‘^(u,y*) := sup R{v,y*]jo,S). (3.3) 

||5||T<e 

In (|3.2I) . Eo{v,y*) is the approximation error and E^’^^{v,y*) is the error gen¬ 
erated by the incompletely known data. The main task is to compute (or esti¬ 
mate) the respective supremum in (|3.3I1 . The difficulty of this task depends on 
the problem and the set of possible data perturbations. 

The described error decomposition is schematically illustrated for a case 
where Eo dominates and vice versa in Fig. [TJ In Fig. [I](top), the approxi¬ 
mation (v,y*) is so far from the exact mean solution {uo,pt) that all admissible 
error measures yield relatively similar value. In this case, it would be reasonable 
to improve the approximation, e.g., by refining the mesh. In Fig. [T] (bottom), 
{v,y*) is so close to {uo,pl) (and the entire set 5(T)) that different admissible 
duality gap functionals provide very different values. This means that the un¬ 
certainty in the data dominates and reducing the error w.r.t. one error measure 
could result in increasing it in another. Since all error measures are equally 
possible, improving the approximation is impossible due to the accuracy limit. 

For numerical analyst, the decomposition (13.21) gives means to compare errors 
generated by different sources. The comparison between different error sources 
can be done, e.g., by computing ratios 


P 


approx __ 


Eo 

^worst 


and 
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P 


indet __ 


^ind 

^worst 


(3.4) 






Figure 1: Error decomposition for a coarse (top) and an accurate approximation 
(bottom). 


Remark 3.2. In simulation practise, is typically of lesser importance 

than . For the sake of completeness, we apply the decomposition /TOP to 

j^hest obtain (analogously to CLTld J 

where 

E^^’^{v,y*) := inf R{v, y*] jo, S) 

\\S\\r<e 

is an error reduction term. It adjusts the value of the approximation error in 
such a way that the value of the error becomes the most favorable one (among 
the possible error measures). Analogously Eo and E™'^, the comparison between 
Eo and also provides information on how close the approximation is to the 
solution set. The approximation is coarse, if Eo ^ and it is accurate 

(close to the solution set), if Eo ~ 

Remark 3.3. Note that is always strictly positive, if the set of admissible 

data is greater than a single value. On the other hand, E^^^^{v,y*) is zero if 
(and only if) {v,y*) belongs to the solution set, i.e., (v,y*) G S(T). In a certain 
sense, and E^^‘^^{v,y*) characterize “the diameter of the solution set” 

diam(„,,.)(5(T)) := F;™-*(n,y*) - F;*’-*(n,y*) 
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measured “from the direction of {v,y*)”. In the special case (v,y*) £ S(T), 
diain(^, resembles more “a radius of the solution set”. 

Remark 3.4. In contrast to Remark \3.S\. it is clear that a real characteristic 
measure of the solution set should be independent of any approximation. It 
should depend only on the problem data and the data indeterminacy. There are 
two natural alternatives for such a measure, 

^worst (y^ ^ p* ) djid sup i?o (<S ( 7 ) ). 

76T 

These two quantities answer to questions “How badly our mean solution may 
violate the necessary conditions related to the family of admissible problems?” 
and “How badly an arbitrary admissible solution may violate the necessary con¬ 
ditions of the mean problem?”, respectively. In linear problems, the duality gap 
is associated with the mixed norm (see Sect. \2.2.1\) and the latter quantity has 
been successfully estimated in m for problems of the form i2.1S\) . 

3.2 Nonlinear reaction diffusion problem 

Consider the problem defined in Sect. 12 . 2.21 where p' = 2,2>p>l, and 
7 ' = 1. Then, V = Hf^(fl), V* = H-^fl), and 

= J dx. 

n 

If the solution is regular enough, then the problem can be written as a stationary 
reaction diffusion problem with a nonlinear reaction term. 

—Au + 7|u|^“^ii =/, inti, (3.5) 

u = 0, on dH. 

Assume that 7 is not fully known, but belongs to a set 

7 GT:={ 7 GL“(r!,R )|7 = 7 o + (5}, (3.6) 

where 70 G L°°(fI,M+) is the known “mean value” and S G L°°(fI,]R) is the 
perturbation. The condition 

||^||L“(n.R) < £ < min7o(x). (3.7) 

guarantees that all elements of T are strictly positive. The compound function¬ 
als (generating the respective duality gap by (12.71) 1 for the problem are 

Df{v, -divy*) = J ^ ~ ~ 

n 

(3.8) 

and 

The following Theorem provides two-sided computable bounds for the indeter¬ 
minacy related error. 
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Theorem 3.1. Let be as in i3.2il . where 7 is defined by and 

i3. 71) . and {v,y*) G x Q, where 

Q ■= {r^* G L‘i(yt, R'^) I divrj* G L^{n, K)}. 

Then, 

where 

E^{v,y*)\= J max{h{-e,v,y*,e),h(e,v,y*,s)}d-K, ( 3 . 9 ) 

n 

£;ind(^^y*) J uia,x{h{-e,v,y*,-e),h{e,v,y*,-e)}}d:xi, ( 3 . 10 ) 

n 


and 

h{d,v,y*,c) := S + (1 - - divy* + /|«) 

+ (i(g - l)(y„ + c)-«-i| - divy* + /r) . 

Proof. The decomposition is derived by linearizing the (nonlinear) term 
containing 7 in (EHl). More precisely, the (first order) McLaurin series of 
z{6) := (70 + i5)^“^ is of interest, i.e., 

z{5) = z{0) + z'{0)5+^z"{c)S^ 

(70 + = 7o^"‘? + (1 - - ig(l - g)(7o + (3.11) 

where c G (0, d) in the remainder term. Applying (13.111) into (13.81) {Dq does not 
depend on 7 ) leads to 

E,{v,y*) = Eo{v,y*) + R{v,y*;io,6), 


where 


R{v:y*-,lo,S) := j (<5(i|?^r + (l-g)7o'’^|-divy* + /|«) 

n 

- d^d(l - g)(7o +c)~'^~^ \ - divy* + /|'?^ dx. 
The remaining task is to estimate sup R{v, y*; 7o, S) from both sides. 

ll'5||j:,°o(fi,R) <E 

This can be done pointwise. Define 

A(x) := i|z;(x)r + i(l - g) 7 o(x)-«| - divy*(x) + /(x)^ 

and 

i?(x; c) := -i(l - ( 7 )( 7 o(x) + c)-«-i| - divy*(x) + f{x)\T 
Then one can write 

Riv,y*'-,lo,S) = J h{S,v,y*,c)dx = J {S{x)A{x.) + S{x.f B{x; c)) dx. 
o o 
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Note that B > 0 a.e. x in fi, since q > I and 70 + c > 0 by the fact that c € (0, S) 
(or c € ((5,0), if 5 < 0) and (13.71) . Thus the integrand is an upwards opening 
parabola w.r.t. S (at each point) and the maximal values are attained at the 
endpoints of the interval (—e,e). Moreover, the function c —5> ( 70 ( 2 ;) + 
is monotonically decreasing. Thus, the maximum of B{x,c) w.r.t. c can be 
estimated from above by setting c = —e and from below by setting c = e. These 
observations lead to (13.9|) and p.lOl) . □ 

3.3 Numerical examples 

Consider the incompletely known reaction diffusion problem described in Sect, 
where := (0,1). The corresponding primal energy functional is 

1 

J{v)= j da; 

0 

and the respective dual energy functional is 

1 

I*iyl = -j{\ I - (y*)' + /I’ + W?) 

0 

The incompletely known function 7 is defined by (13.61) . The approximations Vh 
and are computed by minimizing and maximizing primal and dual energies 
(associated with the “mean data” 70 ) over finite dimensional subspaces, respec¬ 
tively. The minimization and maximization are done by Matlab (version 2013a, 
fminunc -function, quasi-Newton method, BFGS-scheme (Zlllllliaiaol). The 
subspaces are generated by the finite element method using standard Lagrange 
elements. 

Remark 3.5. Generally, there is no need to solve both the primal and dual 
problem. If the applied solution method generates only an approximation of u, 
then the approximation of p* can he typically obtained by postprocessing. 

The accuracy of the pair {vh,y^) is measured using the compound functionals 

1 

DF{v,y') = J (^d7okr + d | - (y*)' + f\'^ - {-(y*)' + /)u^ da; 

0 


and 

Dg{v', -y) = ^llf' + 2 /*lli 2 ((o,i))j 
and the respective duality gap. 

Example 3.1. Let p = 1.2, 70 = 1, and / = 1. Approximations Vh and y'^ 
are computed using linear Lagrange elements. The convergence of the solution 
is shown in Fig. The convergence of the duality gap together with the hounds 
for fTheorem \S.l\) with various levels of indeterminacy in 7 (i.e., e value 
in are depicted in Figure O The bounds for can be easily used to 

generate bounds for ratios and (defined in jS.fD ), which are shown 

in Figure^ If the function 7 is known by accuracy e = 0.005, then refining 
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# nodes 


Figure 2: Example 13. II Convergence of the approximations measured in terms 
of compound functionals and the duality gap 


the mesh would be reasonable. However, if e = 0.01, then all computations 
with more than 60 DOFs are dubious. Moreover, if e = 0.05, then already the 
coarsest approximations with 11 DOFs is sufficient. 

In numerical practise, nonlinear problems have to be solved by some iter¬ 
ative method (for suitable optimization methods, see, e.g., [HIS]). Thus the 
approximation error is generated not only by the mesh, but also due to the 
iteration. For analyst, it might be difficult to judge whether one should refine 
the mesh or improve the iterative solver in order to obtain more accurate ap¬ 
proximations. Moreover, computational cost of the iteration (e.g., optimization 
method) increases quickly as the degrees of freedom are increased and comput¬ 
ing approximations on unreasonably refined meshes (w.r.t. data indeterminacy) 
should be avoided. The following examples show how bounds from Theorem 13.II 
can be applied to define reasonable mesh size using a method of manufactured 
solutions. 

Example 3.2. Let p = 1.2, ‘jo{x) = l + x, and u{x) =sin(27ra;). Thenp*{x) = 
—u'{x) and the respective f{x) is easy to compute by L‘1.5\} . The pair {IhU, IhP*) 
denotes the interpolated solutions on a given mesh. The convergence (using 
linear and quadratic Lagrange elements) of the duality gap together with the 
bounds for E^^^{IhU,LhP*) with different indeterminacy levels e are depicted in 
Fig. O It shows the efficiency of using quadratic elements and reasonable mesh 
sizes for different indeterminacy levels e in id. 71) . In particular. Fig. \^indicates 
that the use of more than approximately 100 DOFs (with quadratic elements) 
is unreasonable for the current problem setting unless the data is known by 
accuracy much less than e = 0.005. 

Example 3.3. Repeat Example 1/1.ill but set u{x) = sin(47ra:). This results 
in a more oscillating f in comparison with Examvle \d.‘A The convergence (for 
quadratic elements) of the duality gap together with the bounds for E™'^{IhU, IhP*) 
with different indeterminacy levels e are depicted in Fig. 
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Figure 3: ExamDle l3.il Convergence of the approximation error in terms of the 
duality gap and the error due to indeterminacy 
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Figure 4: Example 13.11 Bounds for ratios and in (13.41) 
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Figure 5: Example 13.21 Convergence of Eo{IhU, IhP*) for linear (above) and 
quadratic (below) elements together with bounds for E™'^{IhU, IhP*) with vari¬ 
ous indeterminacy levels e in (Id.711 
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Figure 6 : Example 13.31 Convergence of Eo{IhU,IhP*) for quadratic elements 
together with bounds for IhP*) with various indeterminacy levels e in 

(ITTD 


The comparison between Examples 13.21 and 13.31 illustrates that for a more 
complicated problem (highly oscillating righthand side) the reasonable mesh re¬ 
finement level (in comparison with the data uncertainty) is much higher. In 
other words, for a more complicated problem, the accuracy limit dictated by 
the incompletely known data is not met without additional computational ex¬ 
penditure. Note that this procedure could be done in a systematic way by 
precomputing proper mesh refinement levels for a wide range of different known 
u (and the respective /) and indeterminacy levels. Then an analyst could com¬ 
pare his or hers given / and data accuracy to the benchmark computations and 
define a reasonable characteristic mesh size. The same information could be 
used as a stopping criteria for any adaptive solver, since reasonable value for 
Eo{v,y*) would be known a priori. 

It would be possible to study the model problem in great detail, e.g., the 
effect that the parameter p has to the indeterminacy error, compare the inde¬ 
terminacy in the 7 ' and 7 in (I2.14I1 . etc. However, it is not the aim of this 
paper to concentrate on properties of this particular model problem, but to 
demonstrate a convincing methodology that can be applied to study the effects 
of incompletely known data for a wide range of convex minimization problems. 
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